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Analysis of the experimental data on cosmic ray spectra in the framework of the proposed model with two 
types of sources leads to conclusion, that sources with particle generation spectral exponent p ~ 2.85 give the 
major contribution to the all-particle spectrum in the energy range 10 s — 10 7 GeV. 'Fine structure' of spectrum 
around the 'knee' may arise due to presence of nearby supernova type source, accelerating particles up to the 
energies ~ 3 • W 4 Z GeV, if the energy output of such source is ~ 2 • 10 48 erg/source. 



1. Introduction 

Behaviour of primary cosmic ray nuclei spectra 
in the vicinity of the 'knee' is, undoubtedly, one 
of the key elements for solution of problem of high 
energy cosmic rays sources and for establishment 
of acceleration mechanisms. So, for example, 
non-monotonous behaviour of mass composition 
and additional breaks in all-particle spectrum in 
the energy range 10 5 — 10 7 GeV would support 
scenario with supernovae as the major cosmic ray 
sources and acceleration on shock- waves fronts up 
to the energies E max ~ 10 5 Z GeV. Particle spec- 
trum produced by this mechanism may be repre- 
sented as Ssn ~ E~ 2 <d(E max — E), where Heav- 
iside step function 9(x) qualitatively reflects a 
sharp cut-off of the spectrum for E > E max [1,2]. 

The recent results on the proton spectrum of 
Tibet collaboration [3] with no indications on cut- 
off in the spectrum up to ~ 10 7 GeV allow to 
affirm, that proton spectrum in the considered 
energy range is formed by sources with generation 
spectrum being different from Ssn- On the other 
hand, in papers [4,5] it is shown, that all-particle 
spectrum in this energy range has several breaks. 

Reconciliation of experimental results [3] and 
deductions [4] , based on the analysis of E AS size 
spectra and data on Cherenkov radiation of show- 
ers, is possible, in our opinion, only within the 
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framework of the following scenario: 1) particle 
generation spectrum for energies 10 5 — 10 7 GeV 
in the principal cosmic ray sources is differ- 
ent from that of Ssn; 2) contribution of young 
nearby supernova determines 'fine structure' of 
the spectrum and non-monotonous character of 
mass composition behavior in the 'knee' region. 

In this paper we present new calculations of 
energy spectra and mass composition of cosmic 
rays, performed for the above formulated scenario 
of particles generation in two distinct types of 
sources under assumption of anomalous diffusion 
of particles in inhomogeneous galactic medium. 
The main purpose of the calculations is to deter- 
mine conditions under which the proposed model 
(two types of sources with different acceleration 
mechanisms + anomalous diffusion) reproduces 
results of [4,6]. 

2. Diffusion model 

Following to [7,8], cosmic ray transport in frac- 
tal interstellar galactic medium with 'traps' and 
without energy losses and nuclear interactions is 
given by anomalous diffusion equation 

^ = -D(R, a, /3)D^(-A) Q / 2 N(r, t, R) + 

+S(r,t,R). (1) 

Here D = DoR 5 - anomalous diffusion coeffi- 
cient, R — particle rigidity, Dq + — the Riemann- 
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Figure 1. Cosmic ray spectrum in the 'knee' re- 
gion. Dotted lines — data from [4], solid line 
— calculations in the anomalous diffusion model 
with one type of sources S ~ E~ 2 - S5 . 



Liouville fractional derivative [9]: 



astrophysical objects, the solution of equation (1) 
has the form 



N(r,t,R) = 



SoR-p 



T -3/3/a x 



max[0, t-T] 



where function ^3 ^(r) — density of fractional 
stable distribution [10]. 

Analysis of the obtained solution and its com- 
parison with the available experimental data for 
energies below and above the 'knee' allowed to 
determine the basic model parameters [11]. It 
is shown, that for p as 2.85, S as 0.27, anoma- 
lous diffusion model satisfactory describes the ba- 
sic observable features of cosmic ray spectrum 
and mass composition in the wide energy range 
10 2 - 10 10 (see, e.g., [8,11]). The obtained all- 
particle spectrum for sources with S(E) ~ £;-2.85 
practically coincides with 'background' spectrum 
from paper [4] (Figure At the same time 
the calculations demonstrated that for purely 
power-like spectrum S ~ £?- 2 - 85 the model is un- 
able to reproduce additional breaks in the region 
10 5 - 10 7 GeV. 



(— A)"/ 2 — fractional Laplacian ('Riss' opera- 
tor) [9]: 



4,i(a) J \y\ m+a 



where I > a, x € R m , y G R m , 



^/(») = E(- l )*( t )/(*-*v) 



fe=0 



d m ,i(a)= / {l-jv) l \y\- m - a dy. 



If a = 2 and (3 — 1, the equation (1) is just the 
normal diffusion equation. 

For punctual impulse source with power energy 
spectrum S(r,t,R) = S a R- p 6(r)e{T - t)Q(t), 
corresponding to particle generation processes in 



3. Results and conclusions 

To find conditions when the model, formu- 
lated in Introduction, is able to reproduce results 
of [4,6], spectra of nuclei, all-particle spectrum 
and (In A) were calculated with the use of (2) and 
particle generation spectrum in the source of the 
form S(E) = S E- 2 - 85 + SiE~ 2 Q(E max - E). As 
a criterion for determination of particle maximum 
energy, accelerated by supernova type source, we 
used conclusion on the considerable increase of 
(In A) in the region 10 7 - 10 s GeV, obtained 
in [6]. Results of our calculations, given in Fig- 
ures and in Table 1, show, that inclusion into 
system of sources with S(E) ~ £;- 2 -85 adopted 
in [11], of the additional supernova- type source 
(r as 100 pc, t as 10 5 yr), accelerating parti- 
cles up to E nmx as 3 • 10 4 Z GeV with output to 
proton component ~ 2 • 10 48 erg/source, allows 
to describe complicated structure of cosmic ray 
spectrum and mass composition in the vicinity of 
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Figure 2. Comparison of nuclei spectra obtained in the model with two types of sources with experimental 
data (experimental data from [11]). 



the 'knee'. As mentioned above, the all-particle 
spectrum for sources with S(E) ~ iJ~ 2 - 85 practi- 
cally coincides with 'background' spectrum from 
paper [4]. 

In conclusion, also, let us note, that verifica- 
tion of the hypothesis [4] about contribution of 
nearby supernova to the presently observed in 
the Solar system cosmic ray spectrum may be 
accomplished via execution of additional mea- 
surements (or complementary analysis), e.g., of 
H, He and CNO nuclei spectra in some energy 
regions. As follows from our calculations, pre- 
sented in Figure |2 nearby supernova could give 
rise to inhomogeneities in the spectra of these nu- 
clei in the regions of 3 • 10 4 GeV/nucleon (H), 
6 • 10 4 GeV/nucleon (He), 3 • 10 5 GeV/nucleon 
(CNO). Indications on the presence of such in- 
homogeneities in the spectra of He and CNO in 
these energy regions are present in the data of 
JACEE collaboration [13]. 



Table 1 

Mean logarithmic cosmic ray mass (In A) vs. pri- 
mary particle energy in the framework of the 
anomalous diffusion model with two types of 
sources. 

E, GcV (In A) E, GcV (In A) 



10 1 


0.32 


6.3-10 a 


2.05 


2.5-10 1 


0.49 




10 6 


2.13 


6.3-10 1 


0.74 


2.5 


10 6 


2.33 


10 2 


0.90 
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Figure 3. Cosmic ray spectrum in the 'knee' re- 
gion. Comparison of the data [4] with the present 
paper results for the model with two types of 
sources. Dotted line - 'background' spectrum [4], 
dashed line - spectrum [4] with account for young 
nearby supernova, solid line - all-particle spec- 
trum, obtained in the present paper within the 
model with two types of sources. 
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Figure 4. Mean logarithmic cosmic ray mass 
vs. primary particle energy. Dashed areas — ex- 
perimental data from [6, 12], solid line — our cal- 
culation in the framework of the anomalous dif- 
fusion model with two types of sources. 
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